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We study the quantum dynamics of a time reparametrization invariant system with a vanishing 
Hamiltonian. The evolution of the physical degrees of freedom of the system is described, both at 
the classical and at the quantum level, in relational terms by the construction of an internal time 
parameter. We use the Pegg-Barnett phase operator formalism in finite dimensional Hilbert space 
as an essential ingredient. 
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O . I. INTRODUCTION 

:z; 

' The description and the physical interpretation of the quantum dynamics of general covariant systems appear 
[ problematic due to the absence of a clear and unambiguous notion of time for such systems. This lack of a single 
preferred independent external time variable contrasts with the standard formulation of quantum mechanics where 
time plays the role of an absolute element, described by a classical parameter and not by an operator and is required 
for the description of the quantum evolution (it explicitly appears in the Schrodinger equation)^. This tension between 
generally covariant systems and standard quantum mechanics poses both technical and conceptual problems. This kind 
^ l of difficulties are particularly evident when one tries to quantize general relativity. The background independence of 
• the gravitational field, consequence of the fundamental symmetry of the theory, namely general covariance, represents 
^ ] a major obstacle to the application of the standard background dependent methods of quantization of classical fields. 
O^. However techniques have been developed to preserve background independence at the quantum level: Loop Quantum 
Gravity, a candidate for the theory of quantum gravity, realizes a (canonical) quantization of the gravitational field 
in a background independent way Concerning the conceptual issues, several solutions have been proposed in 
order to provide a consistent physical interpretation of quantum generally covariant systems. In this letter we want 
to study one of these proposals advanced by RoveUi [l|, [3, S 0], consisting in describing evolution of the physical 

■ degrees of freedom of the quantum system in a relational manner^. We implement such a relational description in 
a simple Hamiltonian model that shares with generally covariant system the property of being invariant under time 
reparametrization. 

The paper is organized as follows. In the section II we introduce the model of two harmonic oscillators with an 
Hamiltonian constraint and study their properties. In particular we provide a definition of an internal time that allows 
us to describe the classical evolution in a relational manner. In section III we shortly present the canonical quantization 
of the model and its physical Hilbert space (all the details can be found in the references cited) . In section IV we 
translate at the quantum level the classical relational dynamics presented in section II. A fundamental ingredient will 
be the Pegg-Barnett formalism that gives a consistent and well-defined exponential phase operator acting on finite 

■ dimensional Hilbert space. By means of this exponential phase operator we construct well-defined gauge invariant 
d ' operators that represent physical observables of the quantum system and we implement at the quantum level the 

relational evolution of section II. As a by-product of the study of the spectrum of these gauge invariant operators we 
obtain some new relations between the zeros of the Hermite polynomials. 
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^ In the context of field theory, the formulation of the quantum theory uses not only a classical time variable but a whole fixed, non- 
dynamical background spacetime metric, the Minkowski metric. 
^ This proposal is also known as "evolving constants of motion". 
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II. THE MODEL 



With the aim of study the jjroblem of time in quantum general covariant systems, and in particular in quantum 
gravity, Rovelli presented in [l] a simple system consisting of two harmonic oscillators with a total fixed energy^. 
The system has one physical degree of freedom and one first-class constraint The extended phase space T is 
coordinatized by the canonical pairs {qi,Pi);i — 1,2. The symplectic two-form is w = dpi A dqi. The pair (r,w) 
forms a symplectic space. The extended configuration space is defined by the coordinates qt, also called partial 
observables to distinguish them from the complete observables represented by gauge invariant correlations between 
partial observables. We refer to 01 for details. The dynamics, both classical and quantum, may be described using a 
parameter r which is not the physical time, it is not connected with observability but instead it plays the same role of 
the time coordinate in general relativity: physically meaningful quantities are independent from the time coordinate. 
The action of this system has the form 

S[qt,Pi,X]= J dr ^^^pi~ XH{q^,p,)^ , (1) 

where A is a Lagrangian multiplier. This action is invariant under reparametrizations of the parameter r. The 
equations of motion are obtained by the variation of the action with respect to the canonical coordinates {qi,pi), 

dqi _ dH{q^,p^) 

d7 ^ dp, ' 
dpi _ dH{q„pi) 
d7 dq, ' 

while the variation of the action with respect to the Lagrangian multiplier A gives the constraint equation 

H{q,,p,)=0, i = l,2. (4) 

The constraint has the following form 

H{q,,p,) = Hi+H2-M = 0, (5) 

where Hi = ^p\ + \qi is the Hamiltonian of the oscillator 1, and respectively H2 for the oscillator 2, and M is a 
constant. The dynamics is then described by the vanishing Hamiltonian ^ and therefore the system shares with 
generally covariant system like general relativity the important property of being invariant under time reparametriza- 
tions. The physical meaning of such property is that the time parameter appearing in the solution of the equation 
of motion (given below, equations (|6|7p ). namely the parameter r, does not represent a physical time. It is a pure 
gauge parameter. Therefore it is the unfolding of the gauge transformation generated by the first class constraint ^ 
that defines the evolution of the degrees of freedom of the system, i.e. dynamics is gauge. We can parametrize the 
solution of the equation of motion as 

gi(T) = A, sin(r + i = l,2, (6) 

and the momenta are given by 

pj(r) = A, sin(T + (/),), i=l,2. (7) 

These solutions describe an ellipse of radii Ai and A2 and inclination = 01 — 02 in the extended configuration 

space, the plane (qi, 52), with {Ai, A2) G 0, \/2M and (0i, 02) G [0, 27r]. The values of the constants Ai, A2 are fixed 

by the constraint (0]) to satisfy A\ + = 2M. The physical predictions of this system regard the relation between qi 
and 52, not the dependence of qi and 52 on the unphysical parameter r. Consequently the partial observables qi and 
q2 are not physical observables. Instead the physical observables fo the system are given by the so called complete 
observables, namely gauge invariant correlations between partial observables. 



^ We will follow the treatment presented in Q. See also 01 for an alternative description. 
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The Hamiltonian constraint ([4]) defines a constraint surface in the extended phase space T. The restriction of the 
symplectic two-form uj = dpi A dqi to such constraint surface is degenerate. Indeed the Hamiltonian vector field X 
generated by the constraint and given by 

d d 

^^Pta (8) 

OQi opi 

satisfies uj{X) = 0. The integral lines of X on the constraint surface define the "orbits" of lu, namely the motions. 
The constraint surface can be parametrized by the set of independent coordinates {q,pt), where t coordinatize the 
orbits and {q,p) are canonical variables that coordinatize the physical phase space 7. Of course they satisfy {q,p} = 1 
on 7, and they represent physical observables. The symplectic form on j is uj ~ dp A dq. Therefore the solution of the 
equation of motion l|6|7p can be expressed as 

qi = qt{q,p;t), (9) 
Pi=Pi{q,p;t). (10) 

On the other hand, the canonical variables q and p are functions of the coordinates of the extended phase space, 

q = q{qi,Pt), (ii) 
p = p{qz,Pt), (12) 

where the expHcit dependence is fixed by the solution of the equation of motion. The orbit coordinate t as well can 
be expressed as a function of qi and Pii 

t^t{q,,p,). (13) 

This parameter can be interpreted as an internal time with respect to which we describe the evolution of the system. 
The number of such internal time variables is fixed by the constraint: In our case there is only one constraint, so we 
have one parameter t. An interesting property of the internal time variable is that the expression (fT3| is in general a 
multi-valued function, in contrast to the monotonically increasing time variable of non generally covariant systems. 
As a consequence, the values of the variable t can lead to an interpretation of the evolution of the degrees of freedom 
of the system as forward or backward in time. Finally, the relational description of the dynamics of the system is 
encoded in the combination of the solution of the equations of motion ^ and (fTO|) with the definition of the internal 
time variable lfT3|) . 

We will choose the following canonical variables on the physical phase space 

q{q^,P^) = ^{ql+pl)^Al (14) 

p(qi,Pi) = arctan — — arctan — = A0. (15) 
Pi P2 

It is easy to show that and A(j> are conjugate variables, namely {A^, Acf)} — 1. They represent physical quantities: 
They are indeed independent of the gauge parameter t and have vanishing Poisson brackets with the constraint ^ . 
The variable t is given by 

t = arctan — . (16) 

P2 

As already said this variable plays the role of an internal time for the system: It is the parameter with respect to which 
we describe the evolution of the degree of freedom of the system in a gauge invariant way. We notice that the time 
variable iflGj) is not a physical observale, indeed its Poisson bracket with the Hamiltonian constraint does not vanish. 
In particular it results to be equal to 1: {t,H} = 1. Furthermore, the internal time variable iflGj) is a multivalued 
function of 92,^2 because the arctangent only defines t mod 27r. Of course this particular choice of coordinatizing 
the constraint surface is not at all mandatory. We can perfectly choose different variables, for example the couple 



A(j)) with t = arctan -2^. 

If we plug the choice of variables l|14|15p and the internal time lfT6|) into the solution l|9|10p , we obtain 

gi(A?,A(/),<) = ^sin(i + A0), (17) 
92(^2, A0,t) = Jm-AI sin(t- A</.), (18) 
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and the momenta 

Pi{Al,A(f>,t) ^ •y/A2cos(t + A(/)), (19) 
P2{Al,A(l>,t) = ~Aj cos(i- A0). (20) 

The set of equations (|17ll20p give the relational evolution of the Qi and Pi with respect to the internal time (flGl) for the 
particular choice of physical observales (A^, A0). Indeed, choosen a particular combination of a set of variables of the 
extended phase space F as an internal time parameter (in our case the couple ((j2,P2) defining the internal time l|16p ). 
we can describe the evolution of the remaining set of variables, namely qi andpi, as functions of this combination iflGj) 
of 92,^2 for any physical state (A^, A0) of the system. A different interpretation is viable. The expressions (|17ll20p 
give a one-parameter family of physical observables defined on the physical phase space, where the parameter is the 
particular combination of 52 ; P2 given by (flGll . We can then say that (|17H20p are evolving constants of motion in the 
sense of RoveUi [ll,[3,S,[l. 



III. QUANTUM THEORY 



The canonical quantization of the system is straightforward: Upon quantization the constraint ^ becomes a 
quantum Hamiltonian constraint H, and the (classical) Hamiltonians of the two oscillators become positive definite 
operators Hi and H2 ■ The physical Hilbert space Hphys of the system is the space of states solutions of the Wheeler- 
DeWitt equation that defines the quantum dynamics: H\^) ~ 0. This Hilbert space turns out to be of finite dimension 
precisely because of the existence of the Hamiltonian constraint and the positive definiteness of Hi and i?2- It is 
convenient to take the dimensionality of Hphys equal to the integer 2j + 1 = This finite dimensionality of the 
Hilbert space will be important in the following. 

It is useful to introduce the creation and annihilation operators for the two oscillators"', with the usual commutators 
[a, = [b, = 1, [a, b] = [a, 6^] = 0, and the energy eigenstates |ni, 712) where ni and n2 are the numbers of quanta 
of the oscillators qi and 52 respectively, namely ni and n2 are the eigenvalues of the number operators A^i — a^a and 
N2 — b''b respectively. The quantum Hamiltonian constraint H can be expressed in terms of these operators, 

H = Ni+N2 + l-M. (21) 

Therefore the physical Hilbert space is spanned by the vectors satisfying the equation M = ni +712 + 1. It is convenient 
to introduce the quantum number m = -^{ni — 712) that runs from m — —j to m = j, with j = {M — l)/2. So that 
Hphys is spanned by the {2j + 1) states |m) = \j ~ m.j + to). The representation of the states in the space of 
coordinates is 

V'm(9i,92) = (22%(j + TO)!(j -m)!)"^^^i/j+,„(gi)i?j_™(gr2)e"^, (22) 

where Hn{q) id the n-th Hermite polynomial. The dynamics can be obtained from an orthogonal projection operator 
P from the kinematical Hilbert space K. into the physical Hilbert space H, P : IC H, defined by 




and the propagator K of the model is the integral kernel of P. The properties and the classical Hmit of the propagator 
K have been studied in Q. Moreover, following a proposal for computing the graviton propagator presented by 
Rovelli in and making use of the coherent states for this system derived in the two point function has been 
computed in 

Here we adopt a different point of view. In order to describe the relational evolution at the quantum level, we 
want to construct the quantum version of the classical solutions l|17H20p . To do this we will take advantage of a 
property, mentioned above, of the physical Hilbert space, namely its finite dimensionality. This property results to be 
essential for the expression of the quantum operators corresponding to the canonical conjugate variables (A^, A^) of 
the physical phase space of the theory introduced in the previous section. It turns out that such quantum operators 
can be defined in a consistent way. 



■* We use the letter a and b to indicate the annihilation operators for the oscillators 1 and 2 respectively. We take h= 1 in the rest of the 
paper. 
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IV. RELATION EVOLUTION AT THE QUANTUM LEVEL 
A. Complete observables in the Pegg-Barnett formalism 

In a series of papers [HI, Pegg and Barnett constructed a well-defined Hermitian optical phase operator in a 
finite dimensional Hilbert space. Studying the quantum theory of harmonic oscillator in a finite dimensional Hilbert 
space, they succeeded in defining quantum phase states as well as a phase state operator. We use the Pegg-Barnett 
phase formalism to write the Hermitian operator acting on the (2j-|-l)-dimensional Hphys corresponding to the physical 
observable A6. Moreover this formahsm enables us to construct a unitary exponential phase operator exp(±i(/)) (see 
references (iQ. [ill, [l^ for the details of the formalism). Therefore we can consider the following operators, 

qiit) = (e**(7Vi)i/2e-»^^ + e-'*e*^^(7Vi)i/2) , (24) 

I / .--^,1/0 -iA0 _ iA0 /A> Nl/2 



pi{t) = ^ (e**(^i)^ e-*^'^-e-'*e*^'^(A^i)^/"j , (25) 

where Ni is the number operator of the quantum oscillator 1, namely iVi = a'' a. The operator e*'^'^ has the following 
number state representation [l3, [HI, [l3| 

e'^^ = \-j){-J + 1| + I - J + + 2| + . . . + |j - + (26) 

The parameter t appearing in the above expressions is the same parameter defined in section II: The internal time 
given by (fTBl) . 

The main property of the operators qi{t) andpi(i) is represented by their commutations relations with the quantum 
Hamiltonian constraint H |(2T|) . With the expression (|26l) of the exponential phase operator it is easy to show that 
both qi{t) and pi{t) commute with H, 

[qi{t),H]^[p,{t),H]=0. (27) 

These operators are gauge invariant operators, representing therefore physical observables of the quantum system. 
Introducing the operators A = e*^"^ (Ni Y^^ and — (iVi)^/^ g-»A0 ^j^g above operators take the familiar form 

gi it) = i= (e" At + e-^'A) , p, (t) - (e" - e'^M) . (28) 

A and A^ play the role of annihilation and creation operators acting on the finite dimensional physical Hilbert space. 
Their commutation relation is 

[A,A^ = I-{2j + l)\j){j\, (29) 

where / is the identity operator. The difference with the usual canonical commutation relations is a consequence of 
the finite dimensionality of the Hilbert space, and is at the heart of the Pegg-Barnett formalism. It is straightforward 
to calculate the commutation relation between the operators I|24l25p 

Mt),Mt)] = t{AA^-A^A), 

= z(/-(2j + l)|j)(j|). (30) 

Once agian, the above commutator differs from the usual canonical commutator because of the finite dimensionality 
of the Hilbert space. We can write the equations of motion governing the evolution of these operators with respect to 
the internal time t in Heisenberg form as 

^'^^^^^ = ^[Hi'il{t)], (31) 



dt 
dpiit) 
dt 



= i[H[Aii^l (32) 



with H[ = A'^ A ^ Ni. This Hamiltonian differs from the Hamiltonian Hi of the oscillator 1 by a factor 1/2. This 
difference has consequence in the classical Hmit of the propagator, see section IV. C. In the next section we will study 
the spectrum of these operators that turns out to be finite and discrete. 
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B. Spectrum of the complete observables 

The Heisemberg equations of motion for the gauge invariant operators qi (t) and pi (t) , with solution given by the 
expressions (|24|) and ((25|) respectively, represent the equations of motion of a one dimensional harmonic oscillator 
evolving with the time t. The generator of this time evolution is the Hamiltonian H{. Its eigenvalues are the 
eigenvalues of the number operator of the oscillator 1, namely (j + m). We can write the solution of the equations of 
motion adopting the Schrodinger picture: The kets solution of the quantum dynamics have following form 

\qk{t)) ^Nk + m)!)-i/2if,+„(g,)e-«^/2e'(j"+")V), (33) 

where Nk is a normaHzation factor and Hj+m is the Hermite polynomial of order j + m. In order to truncate the 
series at the value m = j, which corresponds to the maximum allowable eigenvalue of Ni, the argument qt of the 
Hermite polynomial must be a zero of H2j+i. There exist 2j + 1 of such zeros (labeled by the index fcV consequently 
the states l|33p span the (2j + l)-dimensional physical Hilbert space. We recover the same result of and |14l. Il5l|. 
The calculation the normalization factor Ni is straightforward: 



{qk{t)\q,{t)) = \N,\' J2 ^"'rni'lk)e-^K (34) 

therefore 



m=0 ■ 



^^ = ^^^^ME=^™('?fe) . (35) 

\m=0 ■ / 

The 2j + 1 states \qk{t)), {k — 1, ...,2j + 1), given by l(33)) are characterized by different properties: 

(i) They are physical states of the quantum system. They are indeed solution of the Wheeler-DeWitt equation, 

H\qk{t)) = 0, (36) 
with H given by (pT]) . In particular the states \qk{t)) are linear combinations of the energy eigenstates |m). 

(ii) They also satisfy a Schrodinger equation describing their evolution with respect to the internal time parameter 
t generated by the Hamiltonian H'l, 

^l\q,{t))^H[\q,{t)). (37) 

(iii) They are eigenstates of the operator qi (t) defined by (|24l) . Their 2j + 1 eigenvalues qk correspond to the zeros 
of the Hermite polynomial of order 2j + 1 . 

(iv) They undergo a cyclic evolution, in fact \qk{t + 2tt)) — \qk{t)). 

(v) They form an orthonormal basis of the physical Hilbert space, indeed 

imm.it)) = Ar,lvl ^^^>^('^^)^^^-(g') - HMH,,^,{qO (38) 

2^J(2j)!(qfc - qi) 

and qk and qi are zeros of the Hermite polynomial of order 2j + 1. Therefore the above expression is zero for 
k I. Consequently we have 

{qiit)\qkit)) = Ski. (39) 
Notice that the right side is the Kronecker 6 and not the Dirac distribution. 

(vi) The transition amplitude from \qk{t)) to \qi{t')) is 

_ 2i 

{qi{t')\qk{t)) = NkNi J2 7^,HMH^{m)e'^''^+''''^^'e^"'^'-''\ (40) 

m=0 



This quantity may receive a clear physical interpretation. The modulus square of l|40|) represents the probability 
of measure the value qi of the physical observable qi at time t' if we have measure the value qk at time t. 
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The property (v) allows to express the physical states in the basis of the energy eigenstates, |m), in terms of the 
physical states \qk{t)), 

2j + l 

1"^) = X! \lkit)){qkit)\m), 

k=l 

2 

= 5^ |g,(i))iV,(2^+'"(j+m)!)-i/2i7^.^„^(5,)e-*e-^(^+™)*. (41) 
fc=i 

Of course the state |m) is independent of the parameter t: The i-dependence of \qkit)) is exactly canceled by the last 
exponential in the above expression. The scalar product between two physical states |m) and |n), expressed in the 
basis of the states \qk{t)), leads to some new relations between the zeros of the Hermite polynomials reported in the 
appendix. 



C. The classical limit 



In order to study the classical limit, we need the expression of the 2j + 1 eigenvalues qt, roots of the Hermite 
polynomial H2j+i, in the limit j ^ oo (corresponding to the classical limit, see Q). Such expressions have been 
derived in the appendix of pJL]: The roots the Hermite polynomial have the following asymptotic form 



We can see that the eigenvalues qk span an interval centered in the origin of the real line of length 0{^/2j + 1). 
The distance between two successive eigenvalues is 0{l/-^2j + 1). When j tends to infinity the set of eigenvalues 
qk becomes dense in the real axis. However, if we naively take the limit j ^ oo, at finite values of the label fc, 
the asymptotic expressions l|42p and (|43|) vanish and consequently they fail to reproduce the correct spectrum of the 
operator qi . This issue has been anaHzed in [3, and the proposed solution consists in a scaHng of the label k of 
the form 



k = ^M±A,^,^ (44) 

TT 

where the constants q and h are real. It has been shown in that the value of b has no infiuence on the result 
obtained in the continuum limit. Its role is to control the convergence of limit to the correct result^. With the 
prescription l(44| any real number q can now be obtained in the classical Hmit. 

It is interesting to study the classical Hmit of the propagator defined in l(40|) . Recalling the Mehler's formula 

tH^i-)Hn{y) ^ = (1 - .^)-/^exp ( ^-V^^-i-^+/>' \ , (45) 

n=0 ■ ^ ^ 

the evaluation of the limit j ^ oo of the propagator l(40|) leads to the expression 

{qi{t')\qk{t)) = NkNle-'^'-''^/^ X 

' ' " ^ (46) 



2sin(t - t') 



exp 



2sm{t~t') 



-{2qkqi-iql+qf)cos{t-t')) 



We recognize in the second line of l(46|) the propagator of a single harmonic oscillator. The presence of the exponential 
factor e~'*^*~* is a direct consequence of the form of the Hamiltonian governing the evolution of the system, namely 



^ The scaling II44II can be modified by the addition of one more parameter useful for the convergence of the classical limit. We refer to 
for the discussion of this point. 
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H'l — Ni. As noted previously this Hamiltonian is not equal to the Hamiltonian of the oscillator 1 because of the 
absence of a constant term 1/2: Hi ~ Ni + 1/2. This difference is responsible for the appearing of the exponential 

g-»(t-t')/2 gg])_ 

D. The two point function 

With the definition of time given above, we can now define a vacuum state associated with the minimum eigenvalue 
of the hamiltonian H[. Of course, a different choice of the time parameter will lead to a different and inequivalent 
vacuum state. This situation reproduces the well know problem of the definition of a vacuum, and more generally 
of a particle concept, in quantum field theory in curved spacetime. In particular the state \m = j) corresponds to 
the lower eigenvalue of H[ . Having a definition of the vacuum state we can compute the two point function for the 
operator qi defined by the values t and t' of the internal time: 

{j\qiit')qi{m (47) 

We recover the same result (apart from a factor 2) obtained in [§|: In particular the calculation of the two point 
function in ^ has been carried out using the coherent states of the quantum system defined in 0. So, taking the 
coherent state corresponding to the vacuum state \ j) and picked on the two points of the classical trajectory (namely 
the ellipse) defined by the two values t and t' of the internal time parameter, the calculation of the two point function 
based on the procedure applied in gives l(47|) . 

V. DISCUSSION 

As mentioned in section II the model of two harmonic oscillators with dynamics fixed by the hamiltonian constraint 
(|4]) was first introduced in [3| in order to study the problem of time for a time reparametrization invariant system. 
In particular in [l| a gauge invariant physical observable describing the evolution of a coordinate in the extended 
configuration space, the plane ((71,(72), say qi, in terms of the coordinate, (72, was explicitly constructed simply by 
eliminating the arbitrary and non-physical parameter r in the solution of the classical equations of motion (|6|7p , and 
interpreting (72 as a real parameter. This defines a set of physical observables qi{t) interpreted as the value of qi when 
q2 has the value t (for every real t). At the quantum level an operator corresponding to qi{t) has been expHcitly 
constructed and it turns out to be only approximately self-adjoint, in particular only when restricted to a subspace 
of the physical Hilbert space called in []| the Schrodinger regime characterized roughly by states described by a wave 
function peaked around a classical trajectory. Mathematically the non self-adjointness comes from a term contained 
in the expression of qi {t) of the form 

v/M + 2L, -t2, (48) 

where the operator = ^{b'^'b— a^a), and a,b are the creation and annihilation operators of section III. For large 
value of t, l(48|) becomes immaginary. This implies that (71 (i) has complex eigenvalues and that evolution can not be 
unitary outside the Schrodinger regime. 

This problem is completely solved by our choice for the operator qi{t) given by (|24|) . which is a well-defined self- 
adjoint gauge invariant operator for each value of the internal time parameter t defined in iflGj) . Of course this internal 
time t differs from the parameter appearing in l(48|) . The operator qi{t) commutes with the Hamiltonian constraint 
and represents therefore a physical observable of the quantum system. Its spectrum has been computed and results 
to be finite and discrete with all the eigenvalues are real. In particular these eigenvalues coincide with the zeros of 
the Hermite polynomial of order 2j + I. These eigenvalues are physical predictions of the theory. 

Moreover, the operator qi{t) realizes the quantum version of the classical quantity qi{Al, A(j),t) defined in (fT7|) . 
We can conclude that the physical observable qi (t) implements the relational description of the quantum dynamics 
of the system of two harmonic oscillator, i.e. expressing the evolution of one oscillator in terms of the other. A key 
ingredient has been the construction of an internal time parameter. Such parameter, given by (fT6|) . depends only on 
the phase space variables of the oscillator 2. It is not a monotonically increasing function as the usual time in quantum 
mechanics. Instead it is a multivalued function of the canonical classical variables ((72,^*2) and as a consequence the 
physical states \qk) undergo a cycHc evolution. Having an internal time, a description of the quantum dynamics in 
terms of a Schrodinger equation satisfied by the physical quantum states of the system is then available, see equation 
(|37l) . We then recover the standard formulation of time-dependent quantum mechanics. However this has to be 
intended in a relational way. Indeed in usual quantum mechanics the phase space variables describing oscillator 2, 
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namely 52 and p2, would have been promoted to the status of quantum operators. Such approach has been adopted 
in [3]. Here we have studied the system in a different way. The variables 92 and p2 are not quantized, they instead 
provide a definition of the internal time. The evolution of the remaining phase space variables qi and pi, both at the 
classical and at the quantum level, is expressed with respect to such internal time. 

Finally, we have obtain as a by-product of the study of the spectrum of the gauge invariant operator l(24|) some new 
relations between the zeros of the Hermite polynomials, see formulas l|A3lA4IA5p . 



APPENDIX A: RELATIONS BETWEEN THE ZEROS OF THE HERMITE POLYNOMIALS 

We have seen that the eigenvectors of the operator qi form an orthonormal basis of the (2j + l)-dimensional physical 
Hilbert space Hphys, and the associated eigenvalues qk are the zeros of the Hermite polynomial of order 2j + 1. Using 
these two properties we derive in this appendix some new relations between the zeros of the Hermite polynomials. 
Starting from the orthonormality of the eigenstates of qi we derived the expression of the energy eigenstates |m) in 
the basis of the physical states \qk{t)) in formula (|4T1) . 

2j + l 2 

H = + rn)l)-^^^Hj+„,{qk)e-^e~'^^+"'^K (Al) 

fc=i 

As already noted, the i-dependence of \qkit)) is exactly canceled by the exponential term e"*'^-'''"™^* in the above 
expression, so that the state |m) is independent of the parameter t. The scalar product between two physical states 
|m) and |n), expressed in the basis of the states \qk{t)), takes the following form, 

2j+l 2 2 

(n|m) = J2 \Nk\' {i'+'^U + my.)-'^'H,+„,{qk)e-^ {2^+'\j + n)\)-^'^ H,+Me-'-^ ■ (A2) 

k=l 

The energy eigenstates |m) form an orthonormal basis of the physical Hilbert space. The scalar product l|A2p reduces 
to Kronecker delta: (n|m) ~ 5m, n. Substituting in (|A2p the expression of the normaHzation factor calculated in 
(|35l) we arrive at 

2j+i / ^-^ 1 \ 

= (2"+" m! n\)-^/^HMHn{qk) ^ (^fe) ■ (A3) 

fe=i \;=o ■ / 

Let's consider this formula for some specific values of m and n. Consider first the case to = n = 0. We obtain 

2j + l 

k=i 1^1=0 2H1.1 yi^^) 



For TO 7^ and n = 0, 



y / , ^'"^'^^^ = 1, TOe{l,...,2j}. (A5) 



k 

The properties of the equations (|A3|A4|A5P will be studied elsewhere 
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